11/0 9/2004 15:05 FAX 6123051228 MUETING RAASCH GEBHARDT @027 

Exhibiti" ' | 

Pavlidis et aL, A Near^tofrared Rision Scheme for Automatic Detection of Vehicle Passengers," 

Proceedings 1999 IEEE Workshop on Computer Vision Beyond the Visible Spectrum: Methods I 

and Applications, 41*48, Fort Collins, CO (June 22, 1999) (cited by the Examiner as XP- 

001028284) . 

A Nesu>Infrared Fusion Scheme for Automatic Detection of Vehicle Passengers 



L Pavlidis 1 * R Symosek 1 B. Fritz 1 

Honeywell Technology Center 

3660 Technology Dr. 

MN65-2500 

Minneapolis* MN 55418 

(pavlidis, symosek, fritz)@htc.honeywell.com 

Fax 612 951-7438 

Tel 612 951-7338 



Abstract 

We undertook a study to determine if the automatic de- 
tection and counting of vehicle passengers is feasible. An 
automated passenger counting system would greatly facil- 
itate the operation of freeway lanes reserved for car-pools 
(MOV lanes). In the present paper we report our findings 
regarding the appropriate sensor phenomenology and ar- 
rangement for the task. We propose a novel system based on 
fusion of near-infrared imaging signals and we demonstrate 
its adequacy with theoretical and experimental arguments. 



1. Introduction 

There are compelling reasons for the existence of an au- 
tomatic passenger counting system in the HOV lane. In 
particular, such a system will be useful in the following re- 
spects: 

1 . It will facilitate the gathering of statistical data for road 
construction planning. The gathering of usage statis- 
tics in the HOV lane is mandated by the U.S. Federal 
Highway Administration. Currently* the gathering of 
data is performed manually. This is obviously labori- 
ous, inefficient, and prone to error. 

2. It will facilitate law enforcement in the HOV lane. 
Currendy* HOV lane enforcement requires substantial 
commitments of State Highway Patrol personnel and 
equipment. HOV lane enforcement has other costs as 
well. These include the risks of high-speed pursuit in 
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lanes adjacent to stop-and-go traffic and the deteriora- 
tion of traffic flow when tickets are issued during peak 
commute periods. 

3. It will enable the States to offer the option to single 
drivers to use some HOV lanes for a nominal monthly 
fee. 

A complete HOV monitoring system suitable for the 
above applications will consist of a passenger detector and 
a license plate reader. Although, substantial work has been 
reported to the technical literature regarding license plate 
readers* work for automated passenger detectors is still in 
its infancy. There are three major technical challenges m 
the development of an automatic passenger detector: 

1. The imaging signal should provide a clear picture of 
the interior of the can The contrast between the human 
silhouettes and the background should be sufficient to 
provide for reliable image processing- 

2. The pattern recognition algorithm that performs the 
passenger detection should exhibit high recognition 
rates and robust behavior. Of course its performance, 
depends to a significant degree on the quality of the 
imaging signal. Even the best pattern recognition algo- 
rithm cannot perform reliably when the imaging signal 
is corrupted with noise, 

3. The system architecture should be designed in such a 
way that will ensure real-time operation, passenger de- 
tection in both the front and back scats, and protection 
from the weather elements. 

In this paper we address the first from the above three 
technical challenges. We describe a novel near-infrared fu- 
sion system thai provides high quality imaging signal both 
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at day and night and in adverse weather conditions. In par- 
ticular, in Section 2 we give an overview and justification of 
our approach. In Section 3 wc describe in detail the theo- 
retical computations that support our assertions. In Section 
4 we present the experimental validation of our hypothe- 
ses. Finally, in Section 5 we conclude the paper and briefly 
mention our ongoing and future work. 

2* Overview 

Our research for a solution to the imaging aspect of an 
HOV system (sensor phenomenology) was guided by the 
following questions; 

1. Is there a band in the EM spectrum that can pene- 
trate through the vehicle's window glass, during day 
and night and, in adverse weather conditions? Do the 
objects of interest (vehicle passengers) have a consis- 
tent appearance in this EM band, irrespectively of their 
physical characteristics? 

2. If there is more than one band* can we fuse the multiple 
bands in a meaningful way to increase the detecting 
power and reliability of the system? 

3. Are there appropriate cameras for these bands that 
have the necessary resolution and speed to live up to 
the requirements of the problem? 
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Figure 1: Electro-Magnetic (EM) spectrum. 

Fig. 1 shows the Electro-Magnetic spectrum- We have 
limited our sensor phenomenology investigation into the 
infrared and visible spectrum regions- Nature- constraints 
our choices below the visible spectrum, since, gamma rays. 
X-rays and, ultraviolet radiation are harmful to the human 
body- Therefore, the typically active systems in these ranges 
cannot be employed in the HOV lane. Technology con- 
strains our choices beyond the infrared region, since mil- 
limeter wave and radio wave imaging sensors are very ex- 
pensive, bulky and* with insufficient resolution [3]. Still, 



the visible plus the infrared range is a huge area of the EM 
spectrum and we had to identify narrow bands within this 
area that are appropriate for ihe task. 

We know from experience as humans that Ihe visible 
spectrum has certain disadvantages for the purpose of the 
particular application- A visible range sensor (like the hu- 
man eye) cannot easily see at night unless is aided by an 
artificial illumination source. Employing a visible range 
flashlight to illuminate the passing vehicles is definitely not 
an option since it will distract the drivers with probably fa- 
tal results. Tinted window glass (now common in certain 
vehicle types) prohibits a clear view of the vehicle's interior 
to visible range sensors (see Fig- 3)- Also, visible range 
sensors are incapacitated during foul weather conditions. 
Finally, passengers produce variable patterns in the visible 
range* depending on their physical characteristics, time of 
day* and illumination conditions. This variability makes the 
machine vision task much more difficult 

From the above discussion, it is apparent that only the 
infrared range held promise for a solution to the problem. 
Within the infrared range two bands of particular interest 
are the reflected infrared (0-7 — 3-0 /im) and the thermal in- 
frared (3.0 - 5.0 /xm, 8.0 - 14-0 jot) bands. The reflected 
infrared band on one hand is associated with reflected so- 
lar radiation that contains no information about the thermal 
properties of materials* This radiation is for the most part 
invisible to the human eye- The thermal infrared band on 
the other hand is associated with the thermal properties of 
materials- We soon found that the thermal infrared band 
was difficult to exploit for HOV purposes because vehicle 
glass severely attenuates EM radiation beyond 2.4 pm (see 
Figs. 2 and 3). 

Fortunately, a major portion of the rcflected-infrared 
range, the so-called near-infrared range (0.7 — 2.4 pm), 
appeared very suitable for the application at hand. In par- 
ticular, we found that; 

1. A camera in this range can safely operate in the HOV 
lane both day and night. During nighttime we would 
need a matching near-infrared illumination source to 
enhance the scene. Provided that the spectral signa- 
ture of the illumination source is deep into the near* 
infrared range, the light will be invisible to me human 
eye. Therefore, no danger of driver distraction exists- 

2. A camera in this range can "see through" both the vehi- 
clc's windshield and its side- windows. The trans mit- 
tance of typical vehicle windows in the near-infrared 
spectrum is at least 40% (see Figs. 2 and 3). 

3. A camera in this range can operate in adverse weather 
conditions. For example, it has been established that 
the near- infrared spectrum is particularly good in pen- 
etrating haze. 



42 



PAGE 28/34 1 RCVD AT 1 1/9/2004 4:00:32 PM [Eastern Standard Time] * SVR:USPTO-ff XRF-1/8 ^ DN!S:8729306 * CSID:6123051228 * DURATION (mm-ss):11-18 



11/09/2 004 15:06 FAI 6123051228 



MUETING RAASCH GEBHARDT 



@029 




Figure 2: Transmittance of a typical windshield. 




Figure 3: 
window. 



4. 



Transmittance of a typical tinted side- 



If the near-infrared range is split into two bands around 
the threshold point of 1.4 titru the lower-band (0.7 — 
1.4 fim) and the upper-band (1.4 - 2.4 /Am), then ve- 
hicle occupants will produce consistent signatures in 
the respective imagery* In th c upper-band imagery, 
humans will appear consistently dark irrespectively of 
their physical characteristics and the illumination con- 
ditions. In the lower-band imagery, humans will ap- 
pear comparatively lighter. This is because human 
skin appears to have very high reflectance just before 
1.4 /im but very low reflectance just after 1.4 t±m (see 
Fig. 4) (2]. 

We found that the intriguing phenomenon of the abrupt 
change in the reflectance of human skin around 1.4 Jim is 
due to the water content of the human body- Water absorbs 
heavily near-infrared radiation above 1,4 pm and thus has 
low reflectance in this region- (see Fig. 5). Humans consist 
70% of water and therefore they exhibit spectral behavior 
very similar to water. Interestingly, other inanimate objects 
in the vehicle scene maintain their reflectance levels almost 
unchanged, below and above the threshold point 1.4 /im. 
For example, see Fjg. 6 for the reflectance diagrams of 




Figure 4: Reflectance of dark skin versus light skih- 
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Figure 5; Reflectance of distilled water. 

some fabric materials commonly found in the interior of 
vehicles. This observation provoked the following hue of 
thought; Ideally, everything but the human skin signature 
should appear proportionally the same in the HOV imagery 
from the two bands- Therefore, by subtracting an image 
from the lower-band from its matching co-registered image 
in the upper-band we can produce a fused image where: 

L The silhouettes of the passengers' faces will be rein- 
forced (big difference) and more clearly stand out 

2. The background will dim away (small difference). 

This increased contrast will facilitate a clean-cut threshold- 
ing of the fused image. The thresholded result will be an 
image where only the face blobs of the passengers remain 
and everything else is eliminated. A good classifier will al- 
ways classify fast and accurately such a simple binary pat- 
tern, ensuring the reliable real-time operation of the HOV 
system. 

The near- infrared camera we found most appropriate 
to use for testing our ideas was the Sensors Unlimited 
SU - 320- In terms of spectral response it was less than 
perfect because it didn't cover the entire spectral range we 
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Figure 6: Reflectance of different fabric materials. 



were interested in (0.7 — 2.4 /mi). Instead, it covered the 
subrange 1.1 — 1.4 for the lower-band and 1.4 - 1.7 /xm 
for the upper-band. The question we had to address, given 
out hypotheses and the particular camera model available, 
was if the signal to noise (S/N) ratio and the speed of the 
camera would live up to the task. The complete set of cal- 
culations and the interpretation of their results are described 
in the next section. 

3 Theoretical Computations 

Our hypotheses, as described in the previous section, 
held great promise. Before, however* we could proceed 
with any experimentation we had to determine if given the 
particular SU — 320 Camera specifications: 

1. We would have had an imaging signal with sufficient 
S/N ratio. 

2. The speed of the camera would have been sufficient to 
capture the vehicle passengers moving at an average 
speed of 65 mph (freeway speed). 

As wc stated earlier, we consider two spectral bards, one 
above the 1.4 \tm threshold point and one below it. We 
assume that two SU - 320 cameras would film simulta- 
neously the same scene. One camera should be equipped 
with an upper-band filter and one with a lower-band fil- 
ter. Both cameras should be equipped with a polarizer dur- 
ing day time to reduce solar glare. They should also be 
equipped with a tele-photo lens. Because of the SU — 320 
camera technical characteristics we limit the upped-band in 
the range 1.4 - 1.7 and the lower-band in the range 
1.1 — 1.4 fim. We will demonstrate our SfN computation 
for the lower-band only, since very similar things also ap- 
ply to the upper-band. The computation of the camera speed 
will take place as pan of che S/N computation. 

The first step in a S/N radiometric computation is to 
determine the amount of irradiation that falls upon the ob- 



jects of interest [1] - in our case the vehicle passengers. 
The spectral irradiance of the sun (our illumination source) 
on a clear day at sea level is approximately ijufms = 
0.008 Watt$/an 2 in the 1.1 - 1.4 pm. wave-band. In 
our computation, however, we consider the worst case sce- 
nario of an overcast day. For an overcast day the irradiance 
value is reduced by 10~ 3 giving a irradiance at The vehicle 
of approximately 

-fpucrccLsft — 10 * I sunny 

= 10"** 0.008 

= SfiWatts/cm 2 . (1) 

We assume that the camera is pointed at the vehicle's wind- 
shield* not at a side-window. The transmittance of the wind- 
shield of a common vehicle in the spectral band of interest 
is approximately 0.4. We assume the worst case scenario of 
a dirty window (see Fig. 2). This results in an irradiance on 
the vehicle passengers of 

Jpassmngcr = 0*4 * lovtrcast 
= 0.4 * 8 

= 3-2 pWattsfan\ (2) 

The second step in a radiometric computation is to de- 
termine how much of the incident irradiation on the objects 
of interest is reflected back to the sensor (the SU — 320 
near-infrared camera in our case). The radiance into a hemi- 
sphere, assuming the worst case skin reradiate of 0-4 (see 
Fig. 4), would be 

•ftpajjen^tr = 0.4 * Ip as jen^cr/fl* 
= 0.4 * 3.2/tt 

= 0.4 fiWatts/cm? - stcr. (3) 

This represents the reflected portion of the passenger irra- 
diation. The rest is absorbed by the passenger's body. The 
reflected radiation has to pass through the windshield, the 
camera lens, the band-pass Alter, and the polarizer to reach 
the near-infrared sensor array. As wc did earlier, wc assume 
a 0-4 windshield transmittance in the spectral band of inter- 
est. We also assume a f/2 camera lens (14-32° cone angle) 
with 0-8 transmittance, a polarizer with 0.4 transmittance, 
and a band-pass filter with 0.6 transmittance. Then, the ir- 
radiance at the sensor array of the SU - 320 camera will 
be 

Camera = 0-4 * 0-8 * 0-4 * 0-6 * 7T * 
Rpattengtr * Sin 2 (14,32 d ) 

= 0-4 * 0-8 * 0-4 * 0-6 * tt * 

0.4* sin 2 (14.32°) 
= 0.006 fiWatts/cm 2 . (4) 
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The SU-320 camera has square pixels with a side of 37.5 * 
10~ 4 cm or an area 



it* 



^ = 37.5 *10~ 4 * 37.5 *HT 4 
= 1.40 *KT 5 cm 2 . 



(5) 



Consequently, the radiant power on ihe camera pixel will be 



= 1.4 *10 -5 * 0.006 
= 0.084 * lO" 1 * Watts. 



(6) 



The camera's detectivity, D\ is D* - 10 13 cm - 
JHzfWatis. The Noise Equivalent Power is re- 

lated to detectivity JP\ pixel area A. and electronic band- 
width A/ by the following equation: 

NBP = (A**f) x '*/D\ (7) 

We already know the value of D* and A- In order to com- 
pute the NEP wc need to also know the value of A/. The 
bandwidth A/ is determined by the exposure time (speed) 
of the camera. In turn, the exposure time depends on the ve- 
hicle speed (v)» the camera's Instantaneous Field Of View 
[IFOV) and range (r), and the footprint of the horizontal 
translation Based on the parameters t*, IFOV, r. 

and f h t wc compute the required exposure lime (speed) of 
the camera such that the image smear is less than 1 pixel. 
Then, we check if the exposure time value falls within the 
operational range of the SU - 320 camera. If it docs, the 
SU - 320 camera is adequate for the HOV task in terms 
of speed. We can substitute the corresponding value for the 
bandwidth A/ in Eq, (7) and continue the process of com- 
puting the SfN ratio. 




from the oncoming traffic. This arrangement ensures that 
the camera is located in a safe place and has the appropri- 
ate field of view. We assume that the camera focuses at the 
centerline of the incoming vehicle, at the level of the pas- 
sengers^ faces OVound «■ 1-2 m). The half width of a 
standard freeway lane is uw« - 1-8 m. We assume that 
the vehicle travels in the middle of the freeway lane. There- 
fore, the lateral distance of the vehicle's centerline from the 
camera is: 

CHMd, - ty„ + uw = 7.5 + 1.S = 8.8 m. (8) 
Finally, wc assume that for a typical vehicle's windshield 
the average width and height are WW = 1 - 5 m 311(5 38 
0.9 m respectively. 

The Instantaneous Field Of View (IFOV) is the camera's 
Field Of View with respect to a single pixel (sec Fig- 8). 
We assume that the distance r is approximately equal to 
ihe camera's range r (r « r — 40 m). Then, the IFOV 
can be computed from the following equation: 

arctan [(/i«,i A /2)/r ] 
arctan [(0,9/2)/40] 



Figure?: Configuration of the camera set. 

Fig. 7 shows the configuration of the SU - 320 cam- 
era relatively to the oncoming traffic. The camera is located 
Cgr<rund — 3,6 m above the ground, Cftttway ~ 7 -5m off 
the edge of the freeway, and at a distance of r — Ark ~~ 



40 m 



IFOV = 



(240/2) 
= 0,0001 rad, 



(9) 



where Jifpa = 240 pixels is the vertical dimension of the 
SU - 320 Focal Plane Array (FPA). 



^ 




Figure 8; Geometry for the computation of a single pixel's 
horizontal translation - top view. 

At time t the camera's IFOV sees a small portion of the 
passenger's face of diameter D. This small face area is what 
is imaged into a single pixel. We can determine the value of 
D from the following equation; 



as IFOV*r 
r= 0.0001 * 40 
= 0.004. 



(10) 
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The angle 6 in Fig. 8 is the angle between the horizontal 
level and the optical axis of the camera. Because we have 
assumed that the camera is focused at the level of the pas- 
sengers's faces (see the geometry in Fig- 7), the angle 6 
is: 

= 3.43° (11) 

Since we know the values for D and 0, from the geometry of 
Fig. 8 we can compute the footprint fht of a single pixel's 
horizontal translation 

fhi = D/sin(9) 

= 0-004/ sm{3-43°) 

= 0-067 m. (12) 

We assume that the passengers travel at the nominal freeway 
speed of v = 65 mph or v = 29.3 m/see. At this freeway 
Speed the footprint fht >s covered at time t / 

tf = fht/V 

= 0.067/29.3 

= 2.28 msec. (13) 

Therefore, the exposure time tcxporure of the camera should 
be t cxparure < 2.28 rneec if we would like to have im- 
age smear of not more than 1 pixel. The operational range 
of the $U — 320 camera in terms of exposure lime is 
127 fjszc — 16-3 msec- Therefore, the required exposure 
time of t €xp oBure < 2,28 msec is within the camera's oper- 
ational range or in other words the speed of the SU — 320 is 
up to the HOV task. We choose to set the exposure time of 
the camera to 1 msec (t tzp04ure = 1 msec < 2.28 msec) 
which corresponds to a bandwidth of A/ - 1 kHz. 

Now, that we have addressed the speed question and we 
know the value of A/ we can substitute the values for j4, 
A/, and D m in Eq. (7) and calculate the NEP 

NEP = 1.18 * 10" 13 Watts, (14) 

Therefore, the camera signal to noise ratio S/N will be ' 

SfN = P vl9a /NEP = 0.7. (15) 

In conclusion, assuming a worst case scenario (overcast 
day, dirty windshield, dark passenger skin) we determined 
that the SU - 320 camera, equipped with a f/2 lens, a 
1-1 — 1.4 iim Alter, and a polarizer, if it is positioned at a 
range of r — 40m from the incoming vehicle and at a height 
of c srovn d = 3-6 m above the ground), will achieve: 



1. An acceptable smear of less than one pixel because 
the required exposure time of 4.44 msec is within the 
camera's speed capabilities. 

2. A poor signal to noise ratio S/N = 0.7. To boost the 
S/N ratio to a higher value in overcast days we need 
to employ an illumination source. This illumination 
source will also be useful during nighttime. If we op- 
erated in the visible spectrum the use of illuminator in 
the HOV lane would be prohibitive, Fortunately, in our 
case the spectral signature of the illuminator should 
match the range 1.1-1.7 /*rn- Since this range is deep 
into the near-infrared spectrum there is no danger of 
distracting the driver and an illuminator for this range 
can be safely employed in the HOV lane. Contrary to 
the overcast sky scenario, the S/N ratio is exception- 
ally good in the case of a sunny day in both bands (see 
Table 1). Therefore, in bright sunlight conditions, the 
system can operate without the requirement for an ar- 
tificial light source. 

4 Experimental Validation 

The above theoretical scenario was subjected to an em- 
pirical test at one of the Department of Transportation traffic 
monitoring and research facilities. The experiment lasted 
for a week in late September 1998. One test lane of the 
freeway, one mile long, was released from traffic and given 
to us for exclusive use. We set up the SU - 320 camera 
with all its accessories above the test lane of the freeway. In 
the absence of a permanent installation device we installed 
the sensor suite in the basket of a cherry-picker. We imple- 
mented the experiment exactly as it was specified in Section 
3. The only deviation was that we used one SU — 320 cam- 
era instead of two for budgetary reasons, lb simulate the 
operation of the two-camera system, we equipped the cam- 
era first with the lower-band filter. We filmed some carefully 
arranged scenes. Then, wc changed the lower-band filter 
with the upper-band filter and filmed again the exact same 
scenes. Later in the lab» wc co-registered the imagery from 
the two bands using a warping transformation. This way, 
we simulated to a certain degree of accuracy the operation 
of two co-registered cameras in the lower- and upper-band 
respectively. 





1-1 - 1.4 pm 


1.4- 1.7 pm 


Sunny Bay 


700-0 


500.0 


Overcast Day 


0.7 


0-5 



Tftble 1: S/N ratios for different day conditions and spec- 
tral bands. 
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Wc used two testing cars that made successive passes 
ihrough the field of view of the near-infrared camera. The 
passes were done at speed increments of 10 rnpK rang- 
ing from 10 - 50 tnph. One of the testing cars was repre- 
sentative of the compact category (Mitsubishi Mirage) and 
the other of the luxury category (Oldsmobile Aurora). The 
experiment bad both day and night sessions. During night 
lime we used an artificial near-infrared illumination source 
in the range 1.0 - 2.0 fim that covered ihe illumination 
needs of both the lower- and upper-band camera versions. 




Figure 9: Near-infrared day time results, (a) Image in 
the band 1.1 - 1.4 /*m. (b) Image in the band 1.4 - 
1.7 /Am. 




Figure 1(H Comparative results between the visible 
spectrum and the near-infrared fusion approach. 

Fig. 9 shows the images from a particular scene in the 
upper- and lower- bands. The image in the upper-band (Fig. 
9(b)) looks in general darker than the image in the lower- 
band (Fig. 9(a)) because the energy of the son illumina- 
tion in the upper-band is less than that in the lower-band. 
Proportionally, however, the face of the vehicle occupant 
looks much darker because of the marked difference in re- 
flectance of the human skin features between the two bands. 
Kg. 10(bl) shows the result of the fusion operation among 
the images of Fig. 9, Fig. 10(b2) shows the result of the 



thresholding operation upon the fused image. As is appar- 
ent in the figure, only the face of the vehicle occupant has 
remained plus a few small noise regions. The noise would 
have not been there if we had two cameras that were ac- 
curately co-registered and operated simultaneously. This 
simple experiment demonstrates the enormous potential of 
data fusion from the two near-infrared bands wc were us- 
ing. Comparatively, Fig- 1 0(a2) demonstrates how much 
more extraneous information remains after a thresholding 
operation in the visible range image of Fig. 10(al). 

5 Conclusions and Future Work 

We have described a novel method to provide high qual- 
ity imaging signals to a system that will perform passenger 
detection and counting in the HOV lane. The method calls 
for two co-regisicrcd near-infrared cameras with spectral 
sensitivity above (upper-band) and below (lower-band) the 
1.4 jim threshold point respectively. The quality of the sig- 
nal remains high even during overcast days and nighttime, 
because we can safely illuminate the scene with an eye- 
safe near-infrared illuminator [4]- The near-infrared cam- 
eras can also provide clear imaging signals even in certain 
foul weather situations, like for example in hazy conditions. 

The crown jewel of the method is the fusion of the co- 
registered imaging signals from the lower- and upper- band 
cameras. Because of the abrupt change in ihe reflectance 
for the human skin around 1.4 ptm, the fusion has as a re- 
sult the intensification of the passenger face silhouettes and 
the diminution of the background. This increased contrast 
allows foT perfect segmentation that leaves in the final pro- 
cessed image only the face blobs of the passengers. Evi- 
dently* such a clean-cut binary imagery will ensure the re- 
liable and fast operation of the pattern classifier that will 
perform the passenger detection task in the future. 

Wc are currently developing the passenger detection al- 
gorithm and results of its performance will be reported in 
another forum. We are also deriving the design for a pro- 
totype version of the HOV system that will be permanently 
installed in an actual freeway site The prototype design 
calls for a slightly different arrangement of the camera set 
than the one described in this paper. Since the HOV system 
should count passengers in both the front and back seats, 
the cameras should not face only ihe vehicle's windshield 
but also part of the side windows. The theoretical and ex- 
perimental results described in this paper are still valid for 
this modified arrangement to a very close approximation. 
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